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1.  INTRODUCTION 


As  the  sophistication  and  complexity  of  modern  digital  communication  systems 
have  increased,  so  has  grown  the  need  for  realistic  and  timely  simulation  of  these 
systems.  This  need  has  provided  the  impetus  for  the  development  of  the  Interactive 
Communications  Simulator  (ICS),  a  flexible,  graphics-oriented,  and  highly  interactive 
hard  ware /soft  ware  system  consisting  of  a  typical  minicomputer  acting  as  host  to  a  fast 
peripheral  array  processor. 

The  ICS  is  continually  under  development  to  reflect  new  trends  in  the 
development  of  communications  algorithms  and  technology.  The  system  has  been  used 
both  to  evaluate  existing  modem  performance  and  to  explore  new  modulation/coding 
and  signal-processing  concepts  pertinent  to  military,  commercial,  and  space  applications. 


The  purpose  of  this  document  is  to  describe  the  models  employed  in  the  ICS  and 
to  provide  a  cursory  tutorial  on  the  general  concepts  of  digital  communications. 


2.  GENERAL  BACKGROUND 


Digital  simulation  provides  a  useful  and  effective  adjunct  to  direct  analytical 
evaluation  of  communication-system  performance.  Indeed,  there  are  situations  in  which 
explicit  performance  evaluation  defies  analysis,  and  results  can  be  obtained  only  through 
either  actual  hardware  evaluation  or  digital  simulation.  The  speed  and  flexibility 
associated  with  digital  simulation  generally  provide  compelling  reasons  for  adopting  this 
approach. 

Unfortunately,  an  accurate  simulation  of  most  communication  systems  of  even 
moderate  complexity  can  be  terribly  time  consuming  on  present-day  general-purpose 
machines.  This  is  due  to  the  large  number  of  repetitive  signal-processing  operations 
that  must  be  performed  in  order  to  obtain  a  statistically  valid  measure  of  system 
performance.  Recently,  a  class  of  fast  floating-point  array  processors  has  become 
available;  this  greatly  facilitates  these  repetitive  signal-processing  operations.  The  use 
of  the  word  “array”  in  this  context  is  intended  to  indicate  a  processor  optimized  for 
handling  large  data  arrays  in  distinction  to  the  large  machines  organized  as  parallel 
arrays  of  independent  processing  elements.  These  array  processors  are  intended  to  be 
employed  as  peripheral  floating-point  processors  in  conjunction  with  a  general-purpose 
host  computer  which  provides  overall  system  control.  In  the  ICS,  the  general-purpose 
host  provides  an  interface  to  the  user  for  graphics  input  and  display  functions.  All 
signal-processing  tasks  are  relegated  to  the  array  processor.  This  hardware 
configuration  provides  a  unique  opportunity  for  accurate  and  statistically  meaningful 
digital  simulation  of  existing  and  future  communication  systems. 

The  ICS  has  exploited  the  potential  of  this  hardware  configuration.  This  has 
resulted  in  an  extensive  hardware/software  system  for  the  digital  simulation  of  arbitrary 
point-to-point  communication  systems.  While  the  system  can  be  extended  to  include 
analog  communications,  the  primary  interest  is  in  digital  communication  systems.  This 
system  makes  extensive  use  of  interactive  computer  graphics  and  includes  a  Digital 
Equipment  Corporation  VAX-11/780  or  PDP-11/40  acting  as  host  to  a  Floating  Point 
Systems,  Inc.,  AP-120B  floating-point  array  processor.  A  block  diagram  of  the  hardware 
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configuration  supporting  the  ICS  software  is  illustrated  in  Figure  2-1. 

The  user  can  configure  a  wide  variety  of  digital  communication  systems  from  basic 
modules  provided  as  system  facilities.  These  facilities  include:  channel  coders/decoders; 
modulators/demodulators  (modems)  for  a  number  of  modulation  strategies;  receiver 
front-end  filtering;  and  various  bit-synchronization  and  phase-tracking  algorithms.  In 
addition  to  additive  white  Gaussian  noise  (AWGN)  channels,  the  simulator  includes  the 
ability  to  model  impulsive  noise  channels  as  well  as  fading-dispersive  channels  typical  of, 
say,  HF  or  tropospheric-scatter  channels.  ICS  also  provides  an  extensive  ECM  signal- 
simulation  facility  as  well  as  an  explicit  diversity  capability. 


2-2 


3.  FUNCTIONAL  DESCRIPTION  OF  THE  ICS 


In  this  section  a  complete  functional  description  of  the  ICS  is  provided.  The 
general  systems  model  is  discussed,  after  which  the  various  usage  modes  of  the  ICS  are 
described.  Finally,  some  typical  graphical  outputs  from  the  ICS  are  discussed. 

3.1  SYSTEM  DESCRIPTION 

The  use  of  an  attached  processor  in  the  configuration  of  Figure  2-1  has  become 
very  popular  for  a  general-purpose  signal-processing  system.  Such  a  configuration  can 
provide  an  excellent  computation/cost  ratio  in  a  variety  of  applications.  The 
availability  of  a  comprehensive  set  of  library  routines  for  such  array  processors  allows 
programmers  versed  only  in  FORTRAN  or  some  other  high-level  language  to  take 
advantage  of  the  peripheral  processor’s  computational  capabilities.  The  use  of  an  array 
processor  in  such  an  environment  involves  three  steps.  First,  data  which  is  to  undergo 
computation  is  transferred  to  the  array  processor.  Second,  parameters  are  passed  to  the 
AP,  along  with  instructions  to  run  the  appropriate  library  function.  Finally,  the 
resultant  data  is  transferred  back  to  the  host.  Depending  on  the  length  of  the  data  set 
transferred,  and  the  number  of  computations  performed  in  the  AP  each  time  this 
process  takes  place,  a  substantial  increase  in  throughput  may  ensue. 

Unfortunately,  the  above  procedure  is  unsuitable  for  the  very  specialized  types  of 
signal-processing  algorithms  involved  in  communications  simulation.  Many  of  the 
functions  are  highly  bit  oriented  and  therefore  unsuitable  for  simulation  using  standard 
library  routines.  Also,  the  overall  computational  efficiency  of  a  system  used  in  this 
mode  is  directly  proportional  to  the  number  of  computations  performed  during  each  run 
of  the  AP  and  inversely  proportional  to  the  amount  of  data  transferred  between  the  two 
computers. 


As  a  result  of  these  considerations,  to  increase  throughput,  all  of  the  specialized 
signal- processing  operations  implemented  in  the  ICS  are  coded  in  the  AP-120B  assembly 
language,  APAL.  This  was  a  decision  made  early  in  the  development  of  the  ICS,  and  it 
has  resulted  in  a  throughput  increase  over  a  comparable  host-based  system  of  about 
three  orders  of  magnitude.  In  the  simulation  mode,  no  data  is  transferred  between  the 
host  and  AP.  Instead,  all  data  is  generated,  processed,  and  analyzed  in  the  AP.  Once 
programs  and  data  have  been  loaded  in  the  initializations  mode,  the  host  supplies  only 
the  minimal  amount  of  control  necessary  to  monitor  the  progress  of  the  simulation.  The 
data  transfer,  processing,  and  control  operations  for  the  initialization  and  simulation 
phases  are  illustrated  diagrammatically  in  Figures  3-l(a)  and  3-l(b). 

The  internal  processing  in  the  AP-120B  has  been  improved  as  part  of  the  ICS 
enhancements  with  the  design  of  an  internal  AP-120B  executive  program.  This  set  of 
routines  has  greatly  simplified  the  interaction  of  modules  and  submodules  within  the 
AP-120B.  The  executive  program  facilitates  efficient  transfer  of  data  and  parameters 
between  various  assembly  language  modules  resident  in  the  AP-120B.  The  executive  is 
called  by  running  programs  to  store  or  retrieve  data  operated  upon  by  other  programs 
without  interaction  with  the  host.  This  feature  has  permitted  a  good  deal  of  code 
standardization  and  has  increased  overall  system  throughput. 

ICS  is  a  highly  graphics-oriented  system.  Through  graphics  modularity,  a  high- 
level  graphics  display  executive  interfaces  to  a  variety  of  graphics  terminals  including 
Tektronix,  Ramtek  and  IMLAC  displays.  Device  independence  is  maintained  to  a  large 
degree  by  employing  low-level  translation  libraries  appropriate  for  each  specific  display. 

3.2  SYSTEMS  MODEL 

A  block  diagram  of  a  generic  communication  system  which  provides  a  model  for 
simulation-software  development  is  illustrated  in  Figure  3-2.  The  information  source 
generally  produces  a  discrete-sequence  or  analog  waveform  which  is  to  be  encoded, 
transmitted  over  a  specified  channel,  reconstructed,  and  delivered  to  a  remote 
destination  or  information  sink.  In  Figure  3-2,  the  purpose  of  the  source  encoder  is  to 
encode  the  source  output,  presumably  in  an  efficient  fashion,  into  a  binary  sequence  {a,} 
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for  subsequent  encoding  and  transmission.  The  source  decoder  at  the  remote 
destination  provides  the  inverse  function.  That  is,  it  employs  the  binary  sequence  {at} 
to  reconstruct  an  approximation  to  the  source  output.  The  binary  sequence  {dj  may 
differ  from  the  actual  transmitted  sequence  {a,-}  due  to  channel  errors.  This  has  an 
obvious  effect  on  the  accuracy  with  which  the  source  decoder  can  reconstruct  the  source 
output. 

Although  source  coding/decoding  schemes  could  easily  be  incorporated  into  the 

simulator,  they  were  not,  because  the  primary  purpose  of  the  ICS  is  the  evaluation  of 

digital  communication  transmission.  As  a  result,  it  is  assumed  that  the  combination  of 

source  and  source  encoder  can  be  modeled  as  a  binary  memoryless  source  as  indicated  in 

Figure  3-2.  The  major  criterion  of  system  performance  is  then  the  bit  error  probability 

Pb,  i.e.,  the  probability  that  fl,  a,.  This  quantity  can  be  evaluated  by 

straightforward  Monte  Carlo  simulation.  There  are,  however,  other  modes  of  usage  of 

this  system.  Each  of  the  remaining  elements  in  Figure  3-2  will  be  described  in  detail. 

The  binary  memory  less  source  is  implemented  as  the  modulo-2  sum  of  a  linear  feedback 

47 

shift  register  sequence  of  length  2  -1  and  a  mixed  congruential  random-number 

47 

generator.  It  has  a  period  exceeding  2  -1. 

Channel  Encoder 

The  purpose  of  the  channel  encoder  is  to  accept  the  binary  sequence  {fl,}  at  its 
input  and,  through  the  insertion  of  controlled  redundancy,  to  produce  the  M-ary 
sequence  {*,-}  at  its  output.  Generally,  the  encoded  output  sequence  is  such  that 
J,- £  0,1  ,...,M-1,  although  in  the  binary  case  it  will  be  convenient  to  assume  that 
£,•  =  ±1,  i  =  1,  2,....  Both  block  and  convolutional  channel-encoding  capabilities  are 
provided  in  the  simulator.  In  the  former,  the  binary  sequence  {a,}  is  segmented  into 
blocks  of  size  k  to  which  n-k  redundant  bits  are  added  to  produce  a  rate  R=k/n  code 
measured  in  units  of  information  bits  per  transmitted  channel  symbol.  There  is 
complete  independence  between  blocks.  Convolutional  codes,  on  the  other  hand,  do  not 
impose  a  block  structure  on  the  input-information  sequence;  they  do  produce  a  sliding 
dependency  over  a  span  of  input  symbols.  The  encoder  can  be  implemented  as  a  shift 


register  into  which  information  bits  are  shifted  ‘b’  at  a  time.  After  each  shift,  n  modulo 
two  sums  are  computed  from  the  contents  of  selected  shift  register  stages.  These  n 
modulo  two  sums  then  constitute  the  encoder  output  sequence  associated  with  a 
particular  branch  of  b-input  information  bits.  The  normalized  code  rate  in  this  case  is 
R=b/n.  Note  that  each  input  information  bit  affects*  [k/b]n  encoder  output  symbols 
where  k  is  the  length  of  the  shift  register.  The  quantity  k  is  called  the  constraint  length 
of  the  code.  A  typical  k=3,  R=l/2  encoder  is  shown  in  Figure  3-3.  Tabulations  of 
good  convolutional  code  constructions  are  available  for  various  k  and  R  (cf.  [1],  [2],  [14]) 
values.  A  description  of  block  codes  can  be  found  in  Appendix  A. 

The  existing  block  and/or  convolutional  codes  are  generally  effective  in  correcting 
random  errors  but  seriously  deficient  in  their  ability  to  correct  burst  errors.  In  several 
important  applications,  the  channel  environments  are  known  to  introduce  burst  errors 
due  to  time-correlated  fading,  multipath,  co-channel  interference,  jamming,  etc.  There 
are  few  constructive  coding  approaches  for  channels  with  memory.  Indeed,  the  few 
existing  approaches  are  extremely  sensitive  to  channel-modeling  assumptions  which  are 
often  gross  estimates.  In  the  vast  majority  of  existing  systems,  the  approach  is  to 
employ  appropriate  interleaving  schemes  in  an  effort  to  make  the  channel  appear 
memoryless.  In  particular,  the  serial  channel  symbol  sequence  {jt,}  is  interleaved  or 
scrambled  in  such  a  fashion  that  successive  channel  symbols  are  transmitted  separated 
by  many  channel-signaling  intervals.  This  interleaving  is  generally  quite  effective  in 
reducing  the  effect  of  burst  errors  and  allows  use  of  the  well-developed  classes  of 
random-error-correcting  codes.  The  interleaver  implemented  in  the  ICS  is  of  the 
convolutional  type  with  a  maximum  symbol  separation  of  30.  It  is  a  particular 
implementation  of  a  type  1  optimal  interleaver  described  by  Ramsey  (20). 


The  notation  [xl  means  the  largest  integer  not  exceeding  x. 
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Modulator j  Transmitter 


The  purpose  of  the  modulator/transmitter  is  to  map  the  channel  encoder  output 
sequence  { J,}  (possibly  interleaved)  into  a  waveform  S  ( t )  suitable  for  transmission  over 
the  channel.  Encoder  outputs  are  presented  to  the  modulator/transmitter  every  Tt 
seconds.  Tt  is  called  the  baud  interval,  or  equivalently  fe—  1/  Tg  is  the  baud  rate.  It 
is  convenient  to  normalize  all  frequency  (time)  domain  quantities  to  the  baud  rate 
(interval).  For  example,  the  bandwidths  of  frequency-selective  channel-filtering  elements 
are  all  normalized  to  fr  The  filter  cutoff,  or  3-dB  point,  would  then  be  specified  in 
digital  frequency  in  units  of  cycles  per  baud.  This  allows  the  user  the  ability  to 
configure  and  execute  a  simulation  program  independent  of  the  actual  baud  rate,  which 
in  turn  facilitates  user  interaction  with  the  communications  simulator  and  provides 
simulation  results  in  a  convenient  parametric  form. 

In  simulating  the  channel-signaling  waveform  s  ( t ),  complex  notation  is 
extensively  used.  This  allows  simulation  of  a  bandpass  signal  in  terms  of  its  lowpass 
complex  envelope  and  eliminates  the  need  for  simulating  a  high-frequency  carrier 
component.  Specifically,  it  is  assumed  that  $  (f)  can  be  expressed  as 

s  (t)  =  y/2  Re  {s  (t)^2*^}  ,  (1) 

where  fc  is  an  assumed  known  carrier  component  and  s  (t)  represents  the  complex 
envelope,  which  can  be  represented  in  the  form 

s(t)  =  se(t)-j  s,{t)  ,  (2) 

where  sc  ( t )  and  Sg  (f)  are  real  quantities  representing  the  inphase  and  quadrature  (I/Q) 
components,  respectively.  By  generating  only  the  I/Q  components,  the  simulator  need 
not  be  concerned  with  the  explicit  details  of  the  RF  portions  of  the  system.  In 
particular,  discrete  samples  of  the  I/Q  components  are  generated  at  equally  spaced 
intervals. 


In  the  construction  of  the  ICS,  it  has  proven  sufficient  to  consider  that  s  (t)  can 
be  expressed  as 

*  (0  =  yPr  Si c\,iuAl-'TcT)  +  J  c2,iuJl-iTrT)\^  .  (3) 

where  T  and  0  represent  random  timing  epoch  and  carrier  phase,  respectively.  The 
receiver  will  be  required  to  estimate  and  track  these  initially  unknown  quantities.  Here 
the  sequences  { and  {C2J  are  determined  from  the  channel-encoder  output 
sequence  {ar,},  while  uc  (t)  and  Ut  (t)  are  elementary  baseband-signaling  waveforms.  By 
appropriate  choice  of  the  mappings  {C^,}  and  {C2i},  together  with  the  baseband¬ 
signaling  waveforms  uc{t)  and  ug  (t),  the  expression  for  S  is  sufficiently  general  to 
include  most  known  modulation  formats.  These  include:  coherent  binary  phase-shift 
keying  (BPSK),  differentially  coherent  phase-shift  keying  (DPSK),  quadrature  phase- 
shift  keying  (QPSK),  noncoherent  frequency-shift  keying  (FSK),  and  quadrature 
amplitude-shift  keying  (QASK).  In  the  case  of  the  offset  or  staggered  quadrature 
phase-shift  keying  (OQPSK),  and  the  constant  envelope  formats  such  as  continuous 
phase  frequency-shift  keying  (CPFSK),  and  the  special  case  of  minimum-shift  keying 
(MSK)  (cf.  [3],  [4]),  this  formulation  must  be  modified  slightly.  Specifically,  in  this  case 

*  (0  =  lci, iuc( <~2 iTrT)  +  j  C2tMt-2iTrT)]  e&  ,  (4) 


where  now  the  elementary  baseband  waveforms  are  translated  by  two  baud  intervals  at 
a  time. 

In  the  case  of  coherent  BPSK,  for  example,  C2 ,  =  0,  ug(t)  =  0,  while 
C\  i  —  (assume  =  ±1)  and  Ue  (f)  is  equal  to  the  pulse-like  waveform* 

*  The  pulse  waveform  Uq  (l)  in  Table  3-1  has  value  unity  over  the  interval  (-  Tf,  Tj  and  tero 
elsewhere. 
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«o  (0  =  *  ;  o  <  *  <  t, 

=  0  ;  elsewhere  . 


(5) 


Similarly,  in  the  case  of  noncoherent  binary  frequency-shift  keyed  (BFSK)  modulation  it 
is  easily  shown  that  Cl  ,•  =  COS©,-  -f  jsin©,-  while  C2  ,•  =  Here  {0,-} 

represents  a  sequence  of  independent  and  identically  distributed  (i.i.d.)  phases  uniformly 
distributed  on  [-7T,  7r).  Without  this  latter  assumption,  the  phase  could  be  estimated  at 
the  receiver  on  the  basis  of  past  transmissions,  thereby  violating  our  assumption  of 
noncoherent  reception.  The  corresponding  baseband  signaling  waveforms  are  given  by 

«c  (0  =  COs(Aa/f/2)  Uq  (t)  ,  (8m) 


and 


u9  ( t )  =  sin(Ao;//2)  %  (t)  , 


(6b) 


where  Au;  is  some  integer  multiple  of  2tt/  Tt  in  order  to  ensure  orthogonality  between 
transmitted  tones.  Additional  choices  of  {C*,}  k=l,2,  and  associated  baseband  pulse 
shapes  uc  (t)  and  us  ( t ),  are  provided  in  Table  3-1. 


It  should  be  noted  that,  in  constructing  Table  3-1,  the  constraint  has  been  applied: 

iT, 


/'*  |5(f)|2  dt  =  E,  ',  t=  1,2,..., 


(7) 


so  that  Eg  represents  the  constant-signal-energy-per-channel  signaling  element  or  baud. 
Indeed,  since  the  simulation  is  to  be  performed  digitally,  replace  (7)  by 

Ni 

2  |s(^A7,,)|2  =  Es ;  t  =  1,2,...  ,  («) 


where  N  is  the  number  of  samples  per  baud  and  is  under  user  control. 


Typical  choices 
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TABLE  3-1 

DESCRIPTION  OF  SOME  TYPICAL  MODULATION  CAPABILITIES 
AVAILABLE  IN  INTERACTIVE  COMMUNICATION  SIMULATOR 


MSK  X2,  X2(M  COS  (7Tl/2Tg)  Uq  (J)|SIN(iri/2Ts)Ug(l-Ts]l  I/O  WAVEFORMS 


are  N=8,  16,  or  32.  The  quantity  A  Tt=  TjN  is,  of  course,  the  sampling  period. 
Observe  that  (7)  is  indeed  satisfied  for  each  of  the  entries  in  Table  3-1. 

The  communications  simulator  then  generates  the  lowpass  I/Q  samples 
{sc(kATg)}  and  {sJ^kATg)}.  For  example,  in  the  case  of  coherent  BPSK  described 
above, 


sc(kATg)  =  x{S/2EgJ T,cos9  ;  {i-l)N  <  k  <  iN , 

while 

sJ(kATg)  =  xiS/2Eg /  Tg sin©  ;  (i-l)N  <  k  <  iN 

corresponds  to  the  ith  successive  baud  interval,  i=l,2...  . 

Spread-Spectrum  Modulation 

Spread-Spectrum  Modulation  is  a  modulation  technique  which  spreads  the  output 
spectrum  of  a  transmitted  signal  over  a  bandwidth  much  greater  than  that  required  for 
the  transmission  of  the  underlying  information  signal.  Spread-Spectrum  modulation  is 
used  for  purposes  of  jamming  resistance,  low  probability  of  intercept,  and  multiple 
access  in  modern  digital  communications  systems. 

By  far  the  two  most  popular  types  of  spread-spectrum  modulation  for  digital 
communications  are  direct  sequence  (DS)  and  frequency  hopping  (FH).  The  DS 
technique  causes  changes  in  the  phase  of  the  transmitted  signal  at  a  chip  rate 
fc  =  1  /Tc  which  is  typically  much  faster  than  the  symbol  rate  fg  =  1  /  Tr  The 
factor  TJTC  or  fc/fg  is  usually  referred  to  as  the  “processing  gain”  of  the  system  and  is 
denoted  as  G.  When  normalized  at  the  baud  or  symbol  rate,  fc  becomes  the  processing 
gain  of  the  system.  The  FH  technique  hops  the  center  frequency  of  the  transmitted 
signal  over  G  hop  slots  or  frequencies.  In  this  case  G  is  also  referred  to  as  the 
processing  gain  of  the  system. 


(«») 


(8b) 
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The  spread  spectrum  function  is  implemented  as  a  multiplication  of  the 
transmitter  output  by  a  spreading  function  dt  ( t ).  Despreading  at  the  receiver  is 
accomplished  by  multiplying  by  a  complementary  function  d2  (f).  These  two  functions 
have  the  property  that 

/0  \  (*)<*2  (0  =  1  •  (10) 

For  the  direct  sequence  system, 

4(0  =  <«.) 

j-  0 

and 

4(0  =  T'-t,),  (lib) 

j-0 

where  p  ( j )  represents  a  pseudonoise  sequence  of  +1,  -1.  The  term  p  (j)  is  generated 
by  a  linear  feedback  shift  register  sequence  (LFSR),  the  precise  nature  of  which  is  under 
user  control.  The  spreading  function  thus  applies  a  biphase  modulation  at  the  chip  rate. 
The  direct-sequence  system  is  illustrated  in  Figure  3-4. 

In  the  FH  system, 

4  (0  =  E  V  M  -  tj  -  r„)e-«'  (m) 

7=0 

and 
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(!»>) 


<k(t)=  S  u-i(l  ~  tj  ~  Tu)e  1  ■ 

Here  a )j  is  a  pseudorandom  frequency  which  is  chosen  at  the  start  of  each  hop 
interval  tj  and  lasts  for  the  duration  of  the  hop,  Tit 

The  processing  gain  G  and  the  hop  time  T y  are  under  user  control.  A  block 
diagram  of  the  frequency-hopped  system  is  shown  in  Figure  3-5. 

Because  of  array-processor  memory  limitations  of  32K,  the  spread-spectrum 
processing  gain  is  limited  to  G  =  512. 

Finally,  one  last  function  provided  in  the  modulator/transmitter  module  needs  to 
be  mentioned.  In  particular,  any  narrowband  filtering  or  nonlinear  amplifier  distortion 
effects  will  be  simulated  as  part  of  the  modulator/transmitter.  For  example,  if  h  (<) 
represents  the  impulse  response  of  a  narrowband  or  bandpass  filter,  it  can  be  expressed 
in  the  form 

h(t)  —  2 Re{h  ( t )e/2*fct}  ,  (is) 

where  h(t)  —  hc(t)  -  jhJJ)  represents  the  baseband  equivalent  impulse  response.  Here 
hc(t)  and  hj[t)  are  real  lowpass  functions  representing  inphase  and  quadrature  filtering 
effects,  respectively.  Let  Sq(  t)  represent  the  complex  envelope  at  the  output  of  this  filter, 
excited  by  S(t)  at  its  input.  That  is,  Sg(f)  =  s(t)  *  h(t)  where  *  represents 
convolution.  Expressing  50  (t)  in  terms  of  its  inphase  and  quadrature  components 
according  to 


«0  (0  =  S0c  (0  -  j's0»  (0  > 


(14) 


«0c  (0  =  *e  (0  *  K  (0  -  *,  (0  * h,  (0  , 


(15*) 


and 


«o,  (0  =  *c  (0  *  K  (0  +  s,  (0  *  he  (<) ,  (isb) 

which  is  illustrated  more  conveniently  in  Figure  3-6.  A  variety  of  filtering  options  are 
available  to  the  user  in  the  ICS.  These  filters  are  implemented  as  digital  filters  and 
include  such  options  as  linear-phase  finite  impulse  response  (FIR)  filters  and  a  variety  of 
infinite  impulse  response  (IIR)  filters  such  as  maximally  flat  Butterworth  response  and 
Chebyshev  equiripple  response  filters.  Filter  parameters  such  as  cutoff  frequency 
(normalized  to  the  baud  rate),  roll-off  characteristics,  passband  phase,  and  amplitude 
characteristics  are  chosen  by  the  user  from  an  appropriately  formatted  menu  list.  While 
it  would  be  possible  to  provide  various  nonlinear  saturation  effects  associated  with  the 
RF  power  amplifier  in  the  modulator/transmitter  module,  such  a  capability  is  not 
present  in  the  ICS  at  this  time. 

Channel 

In  choosing  a  channel  model  for  the  ICS,  it  was  important  to  provide  a  fairly 
general  approach  which  would  encompass  the  entire  range  of  potential  applications.  A 
useful  channel  model  felt  to  satisfy  this  requirement  consists  of  a  three-component 
fading-multipath  channel  plus  additive  noise.  In  particular,  the  fading-multipath 
channel  is  assumed  to  consist  of  a  single  direct  path,  a  specular  multipath  component, 
and  a  diffuse  multipath  component,  as  illustrated  in  Figure  3-7.  Here  S  (t)  represents 
the  complex  envelope  of  the  transmitted  signal  component,  which  is  assumed  to  be 
expressible  in  terms  of  its  in  phase  and  quadrature  components  according  to 
s  ( t )  =  sc  ( t )  -  js8  (/).  The  quantities  and  ^  in  Figure  3-7  represent,  respectively, 
the  specular  and  diffuse  signal  energies  normalized  to  the  energy  in  the  direct  path  so 
that  energies  in  the  three  components  sum  to  1.0.  These  parameters  are  under  direct 
user  control.  In  addition,  the  user  is  allowed  to  adjust  the  relative  phase  <f>  of  the 
specular  multipath  component  relative  to  the  direct  path.  Other  parameter  choices 
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include  the  nominal  differential  delay  T0  and  the  differential  doppier  f0  associated  with 
the  two  multipath  components.  The  channel  model  assumes  the  same  Nominal 
differential  delay  and  differential  doppier  for  the  two  multipath  components.  The  delay 
and  doppier  spreads  about  these  nominal  values  will  be  determined  by  choice  of  the 
channel  scattering-function  <r(r,  f)  associated  with  the  diffuse- multipath  component.  In 
this  context,  the  quantity  <j{t,  f)  represents  the  energy  associated  with  delays  in  the 
range  (r,T  -f  dr)  and  doppier  shifts  in  the  range  (/,  /  -+-  df). 

The  diffuse  multipath  component  has  been  implemented  by  means  of  the  tapped 
delay-line  model  illustrated  in  Figure  3-8.  This  is  a  fairly  general  model  for  fading 
dispersive  channels  under  the  so-called  wide-sense  stationary  uncorrelated  scattering 
(WSSUS)  assumption  (cf.  [5],  (6]).  The  delay  T  in  Figure  3-8  is  assumed  to  be  equal  to 
the  Nyquist  interval  7=1/  W,  where  W  is  the  transmitted  signal  bandwidth.  In 
practice,  a  sampling  interval  somewhat  shorter  than  this  may  be  desirable  in  some 
applications. 

The  quantities  gi  ( t )  ,  t=0,  ±  1,±2,  .  .  .  ,  ±  N ,  in  Figure  3*8  represent 
tap-gain  functions  which  are  modeled  as  mutually  independent  complex  zero-mean 
widesense  stationary  Gaussian  processes.  In  the  communications  simulator,  all  the  tap- 
gain  functions  are  restricted  to  possess  second-order  Butterworth  power  spectral 
densities  of  the  form 


B? 


Bt  +  f  ' 


*  =  0,±1, 


,±N, 


(is) 


where  B,  represents  the  3-dB  bandwidth  in  Hz.  Specification  of  the  3-dB  bandwidths 
B,-  string,  then,  effectively  determines  the  doppier  as  a  function  of  delay  profile  as 
illustrated  in  Figure  3-8.  Similarly,  specification  of  the  multiplying  factors 
<7j,  i  =0,  ±1,  ±2,...,  ±/V  determines  the  delay  profile.  The  parameters  N,  B,-,  and 
<7j,  i  =  0,  ±1,  ±2,...,  ±N  are  specified  by  the  user.  The  tap  gain  functions  are  then 
obtained  as  the  output  of  appropriately  defined  digital  filters  excited  by  white  Gaussian 
noise  inputs.  While  other,  more  general  doppler-as-a-function-of-delay  profiles  are  easily 
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incorporated  into  the  ICS,  the  approach  described  here  is  expected  to  prove  adequate  for 
the  intended  applications. 

Note  that  the  channel  model  as  illustrated  in  Figure  3-8  includes  provision  for  an 
additive  noise  component.  An  adequate  model  for  a  wide  range  of  applications, 
including  ELF/VLF  and  HF,  can  be  represented  as  the  linear  combination  of  a  low- 
density  shot  noise  process  and  additive  white  Gaussian  noise  (AWGN).  That  is,  the 
complex  envelope  n  (f)  of  the  additive  bandpass  noise  can  be  expressed  as 

n  =  y  (<)  +  w  (f)  .  (17) 

Here  y  (t)  represents  a  complex  low-density  shot  noise  process  while  w  [t)  represents 
complex  AWGN  with  double-sided  noise  spectral  density  N0/2  watts/Hz.  Low-density- 
means  that  the  interarrival  times  of  the  impulses  are  relatively  long  compared  to  typical 
impulse  durations.  The  low-density  shot  noise  component  y  (t)  accounts  for  the 
relatively  infrequent  high-level  and  time-resolved  impulse  noise  hits  due  to  atmospheric 
discharge  phenomenon  or  various  contributions  due  to  man-made  noise.  The  Gaussian 
noise  component  w  ( t),  on  the  other  hand,  is  due  to  a  combination  of  front-end  noise 
and  the  large  number  of  low-level  and  overlapping  impulse  hits.  The  latter  contribution 
is  expected  to  exhibit  Gaussian  behavior  from  central  limit  theorem  considerations. 

The  low-density  shot  noise  component  can  be  modeled  as  the  output  of  a  linear 
time-invariant  filter  excited  by  an  amplitude-modulated  impulse  train  or  a  point  process 
at  its  input.  The  filter  generating  the  complex  low-density  shot  process  y  (t)  is  assumed 
to  possess  impulse  response  h  ( t ),  or,  equivalently,  system  transfer  function  H  (s).  A 
block  diagram  of  the  channel-noise  model  is  illustrated  in  Figure  3-0.  The  impulse 
response  h  (t)  will  in  general  be  complex,  possessing  an  inphase  component  hc  (f)  and  a 
quadrature  component  hf  ( t ).  As  indicated  in  Figure  3-9,  the  input  to  the  linear  filter 
generating  the  low-density  shot  noise  component  y  (t)  is  U  (<),  assumed  to  be  of  the 
form 
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(18) 


U  (t)  =  . 

*=0 

Here  { ttj  is  a  sequence  of  complex  independent  and  identically  distributed  (i.i.d.) 
random  variables  which  can  be  expressed  in  terms  of  I/Q  components  according  to 

«,  =  ««  -  ;  *=o,i,  •  •  d») 

The  process  (N(t),  f>0}  in  (18)  is  a  point  process,  in  particular  a  counting  process, 
whose  event  times  are  described  by  the  sequence  {f,}.  More  specifically,  (N(t),  <>0} 
undergoes  a  unit  step  at  each  of  the  event  times  {fj.  The  Poisson  process  is  a  good 
example  of  a  counting  process.  The  communications  simulator  possesses  the  ability  to 
simulate  stationary  renewal  counting  processes  possessing  Gamma-distributed 
interarrival  times  with  p.d.f. 


/(*)  =  exp {-*/&}  J  *>0  ,  (20) 

r(u)py 

where  t/>l,  and  the  parameter  /?  is  defined  according  to  /?  =  1  /t/\  with  X>0  fixed. 
The  quantity  X  represents  the  average  number  of  events  per  unit  time.  This  quantity  is 
specified  by  the  user  in  units  normalized  to  the  baud  duration,  Tg  seconds. 
Specification  of  the  two  parameters  U  and  X,  then,  completely  defines  the  point  process 
(N(t),  f>0}. 

For  example,  if  t/=  1,  then 


/(/)  =  Xc  Xz  ;  j>0  .  (2i) 

This  is  the  exponential  distribution.  In  this  case,  {N(t),  t>0)  is  Poisson,  with  X  the 
average  number  of  hits  per  unit  time.  Another  special  case  is  obtained  by  letting  ls—>00 
while  holding  X  fixed.  The  result  is 


3-25 


f 


(22) 


f(x)  =  ^2-1/X) 

where  £(')  is  the  delta  function.  In  this  case,  the  impulses  are  periodic  with  rate  X  per 
unit  time.  By  introduction  of  the  class  of  driving  point  processes  described  here,  a 
parameterized  impulsive  noise  model  is  provided  which  includes  these  two  extremes  as 
special  cases.  The  class  of  point  processes  and  hence  tt  (f),  is  easily 

generated  on  a  digital  machine. 

It  is  assumed  that  the  I/Q  components  { and  {  ttjJ  can  be  described  by 

Uci  =  /?,<X> S0  i  ,  (28a) 

and 

uH  =  Rp in©  t  ,  (23b) 

with  {/?,}  an  i.i.d.  sequence  possessing  common  probability  density  function  (p.d.f.) 
while  m  is  likewise  an  i.i.d.  sequence  uniformly  distributed  over  [  — TT,  7T  ].  Related 
work,  [7]  -  [8],  on  modeling  impulse  noise  phenomenon  indicates  that  an  appropriate 
choice  of  p.d.f.  fX  }  is  the  lognormal  distribution 

m  =  ;  ’  ,2,) 

where  fi  and  a2  are  the  mean  and  variance,  respectively,  of  a  Gaussian  variate  g  for 
which  r  =  e9.  Another  useful  choice  is  the  power-Rayleigh  distribution  defined  by 

MR)  =  ~7^’Ra  lexp{~(R/ ro)a)  ;  R>°  *  f25) 

Rq 

where  the  characteristic  exponent  a  lie3  in  the  range  0<cr<2,  and  R0  b  a  scale 
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parameter  to  be  specified.  The  sequence  (uj  is  thus  easily  generated  digitally.  Some 
guidance  on  choosing  parameter  values  to  match  observed  impulse  noise  characteristics 
is  provided  in  [7j  -  [8].  These  studies  indicate  that  values  of  a  between  0.2  and  0.7  for 
the  power-Ray leigh  case  allow  the  impulse  statistics  to  match  observed  impulse 
phenomena. 

The  low-density  shot  noise  process  y  (f)  appearing  at  the  output  of  the  linear  filter 
in  response  to  u  ( t)  at  its  input  can  be  expressed  as 

y  (0  =  yc  (0  -  iy,  (0  -  (26> 

where  the  I/Q  components  yc  (f)  and  y8  ( t )  can  be  generated  according  to 

yc  {t)  =  Ri  *  W *-*,)coser/i,( M.JsineJ  ,  (27a) 

*=0 


y*  (*)  =  If  Ri  *  [ hc{t-tdcos9i  +  MK)sin0J  .  (m>) 

t=0 

Here  hc  ( t )  and  hg  (f)  are,  as  described  previously,  the  I/Q  components  of  the  complex 
impulse  response  function  h  ( t ).  The  terms  he  (£)  and  ht  ( t )  are  both  assumed  to  be 
equal  to  the  same  impulse  response  functions  h$  (f).  Several  choices  for  (t)  are 
possible.  For  example,  the  simplest  choice  is  the  pulse  waveform  of  duration  tc  seconds 
normalized  to  the  baud  duration  as  illustrated  in  Figure  3-10.  The  corresponding 
impulse  response  is 

M<)  =  «-i  (0-«- lit-tc) . 

where  U  ,  (t)  is  the  unit  step  function 


(28) 


{o ;  <<o 


(29) 


The  duration  tc  is  specified  by  the  user  such  that  fc«l/ lambda,  so  that  the  impulse 
process  is  truly  low  density.  Other  choices  of  pulse  waveform  available  are  illustrated  in 
Table  3-2. 

Finally,  the  AWGN  component  which  is  added  to  the  low-density  shot  noise 
component  to  complete  the  channel  noise  model  is  described  in  terms  of  its  I/Q 
components  according  to 


w{t)  =  we(t)-jw,{t)  . 


(30) 


The  components  wc  ( t )  and  wg  (f)  will  be  modeled  as  mutually  independent  zero-mean 
white  Gaussian  noise  processes,  each  possessing  double-sided  noise  spectral  density 
Aq/2  watts/Hz.  These  components  are  easily  simulated.  Gaussian  variates  are 
accurately  generated  in  the  ICS  by  generating  a  Rayleigh  random  variable  R  and  a 
uniform  variate  0  on  the  interval  [  -7r,7T  ].  The  transformations 
gi=R sin#  and  g2—Rcos$  then  produce  two  independent  Gaussian  variates.  This 
procedure  is  much  more  accurate  than  approximating  the  Gaussian  distribution  with  a 
sum  of  uniform  variates. 

Generic-Channel  Design 

The  ICS  provides  the  ability  to  design  high-frequency  (HF),  tropospheric-scatter 
(TROPO),  or  line-of-sight  (LOS)  generic  transmission  channels  in  terms  of  parameters 
which  are  physically  meaningful  for  each  channel.  A  translation  program  accepts 
generic  link  information  for  each  of  the  three  links  and  generates  the  appropriate 
fundamental  parameters  for  the  ICS  fading/impulsive  channel  models.  Each  of  the 
channel  design  programs  is  described. 
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TABLE  3-2 

PULSE  WAVEFORMS  AVAILABLE  IN  IMPULSE 
NOISE  GENERATOR 


PULSE  TYPE 

IMPULSE  RESPONSE 
hQC0 

SYSTEM 

TRANSFER 

FUNCTION 

H0(s) 

UNIPOLAR  PULSE 

u.,(t)  -u.|  Ct-tc) 

he  *stc 
s 

EXPONENTIAL  PULSE 

e'alu  (0 

1 

s  +  a 

DECAYING  SINE 

1  e~at  SINatctu_i  (0 
c 

1 

(S  *a)2  Hi )c* 

DECAYING  COSINE 

e'a,COS  w  ctu.,CU 

s  +  a 

Cs  +  a)* 

The  LOS  channel  is  usually  modeled  as  a  direct  path  with  additive  noise.  The 
pertinent  parameter  to  be  determined  is  the  received  SNR,  which  is  calculated  as: 

Eg/ N0=  P Tjbu  -  L OSSiB  -  NdBm  -  10 LOGBr  ,  (31) 

where  the  path  loss  LOSSdB  can  be  expressed  as 

LOSSdB  =  -G 

TiB~GRdB  +  36.6  +  20 LOG(FMBz)  +  20 LOG(Dmiles)  .  (32) 


These  equations  depend  on  the  following  parameters: 


Transmitted  Power/dBm 
Path  Distance  (Stature  Miles) 
Receiver  Noise  Floor  (dBm) 
Transmitter  Antenna  Gain  (dB) 
Receiver  Antenna  Gain  (dB) 
Baud  Rate 
Frequency  (MHz) 


-  Pj-dBm 

'  Emiles 
■  NdBm 

-  GTdB 

*  GRdB 

-  Bb  (bits/sec) 

'  E MHz  • 


The  EJN0  thus  calculated  may  then  be  entered  in  either  the  VALIDATION  or 
SIMULATION  modes. 

The  TROPO  interface  program  is  based  on  a  model  developed  by  the  Defense 
Communication  Engineering  Center  (DCEC)  and  is  based  on  the  tapped-delay-line 
model  described  in  [19j.  The  first  pass  of  the  program  computes  a  power-versus-delay 
profile  after  the  user  inputs  parameters  related  to  transmission  distance,  antenna 
beamwidths,  and  transmission  frequency. 

The  resulting  powers  at  discrete  delay  intervals  become  the  tap  variances  for  the 
ICS  fading  dispersive  channel.  The  second  pass  of  the  TROPO  design  routine  calculates 
received  average  Ej N0. 
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In  addition  to  the  parameters  needed  to  compute  the  tap  gain  variances,  the  user  must 
also  input: 


Transmitted  Power  (dBm)  -  PfdBm 

Receiver  Noise  Floor  (dBm)  -  NdBm 

Transmitter  Antenna  Gain  (dB)  -  GTdB 

Receiver  Antenna  Gain  (dB)  -  GBdB 

Frequency  (MHz)  -  E mHz 

F[Bd]  attenuation  function  (dB)  -  F\6d)  see  [3] 


Atmospheric  absorption  factor  [4j  -  Aa  . 


The  attenuation  function  F\6d)  determines  the  path  loss  as  a  function  of  the 
product  of  the  scatter  angle  $  and  great  circle  path  distance  d.  The  atmospheric 
absorption  factor  Aa  is  a  function  of  frequency  and  path  length.  Both  of  these 
quantities  are  usually  read  from  tables  or  graphs. 

EJ N0  is  then  calculated  as  follows: 

FJ N0  —  PrdBm  ~  Fo$sdBm  -  NdBm  -  10 LogBR,  (33) 


where 


Foss  dBm  —  30logF  mhZ  -  2Ologd0m  -I-  F[d)  +  Aa  -  GTGr  (34) 


The  HF  channel  is  typically  modeled  as  a  one-path,  two-path,  or  three-path 
dispersive  channel  with  additive  impulsive  noise.  The  parameter  translation  program 
employed  in  ICS  is  based  on  the  DICEF  HF  MEDIA  SIMULATOR  developed  by 
SIGNATRON,  Inc.  [20).  The  user  specifies  the  following  information: 

-  BR  bits/sec 

-  1,  2,  or  3 

-  D». 


o  Baud  rate 
o  Number  of  paths 
o  Delay  per  path 
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o  RMS  doppler  spread  per  path 
o  Power  per  path 
o  Mean  doppler  shift 

For  the  impulse  noise  (assumed  rectangular)  the  user  specifies: 

o  number  of  impulses  per  minute 
o  width  of  each  pulse 
o  periodic  or  random  arrivals 
o  power  in  impulsive  component. 

Explicit  Diversity 

The  effect  of  fading  is  often  minimized  by  the  use  of  space-diversity  or  frequency 
diversity  techniques  as  illustrated  in  Figure  3-11.  These  techniques  are  based  on  the 
hypothesis  that  simultaneous  fading  on  both  radio  transmission  channels  is  unlikely.  In 
a  frequency-diversity  configuration,  the  same  information  is  fed  to  two  (or  more) 
transmitters,  TX  1  and  TX  2.  A  wide  frequency  separation  ensures  less  correlation 
between  the  fades  of  the  individual  signals;  therefore,  improved  performance  may  be 
realized.  In  the  space-diversity  system,  the  same  frequency  band  is  used  in  transmission; 
however,  the  two  receiver  antennas  are  separated  vertically.  This  ensures  that  the 
signals  travel  through  different  transmission  paths  and  are,  thus,  not  likely  to  be 
affected  by  fading  in  the  same  way.  Combinations  of  both  space  diversity  and 
frequency  diversity  have  also  been  used  to  achieve  increased  performance.  An  explicit 
diversity-simulation  facility  has  been  incorporated  into  the  ICS. 

Two  popular  methods  of  combining  the  outputs  of  diversity  channels,  selection 
combining  and  maximal- ratio  combining,  have  been  implemented.  Each  method  uses 
SNR  information  in  a  different  way  to  combine  the  channel  outputs.  A  signal-to-noise 
estimate  is  derived  as: 


(35) 


SNR  = 


^(n)  ’ 


where  r  (n)  is  the  short-term  average  signal  level  and  <T2  (n)  is  the  noise  variance. 
Both  quantities  are  evaluated  at  sample  time  (n).  The  signal  level  estimate  is  given  by: 


r(n)  =  -j-  E  lr  (•  )l  • 

11  r  i=r»-nr+l 


where  nr  represents  the  number  of  channel  output  samples  r  (*)  that  are  to  be  used  in 
the  estimate.  The  variance  estimate  is  coupled  to  the  signal  estimate  and  is  given  by: 


'2(«)  =  ^~  E  (|r(.)|-r(»))2, 


a  t==n~n„+l 

where  na  is  the  number  of  samples  used  in  the  estimate. 

The  selection  diversity  algorithm  accepts  data  from  an  arbitrary  channel  among 
the  D  diversity  channels  at  the  start  of  transmission  and  switches  to  a  different  diversity 
path  if  the  new  path  has  the  highest  SNR  estimate  of  D  paths  and  if  this  SNR  exceeds 
the  SNR  of  the  current  path  by  an  amount  A. 

Maximal-ratio  combining  accepts  data  from  each  of  the  D  channels  and  weights 
the  data  of  each  channel  by  its  corresponding  SNR  estimate  before  adding  to  produce  a 
final  output  r1  (i): 


’  (0  =  E  rM/SNRAi). 


Electronic  Countermeasure  (ECM)  Simulation 


A  variety  of  Electronic  Countermeasure  Simulation  modules  have  been  included  in 
the  ICS.  The  jamming  waveform  capabilities  include  multitone-modulated  AM,  noise- 
modulated  AM,  multitone-modulated  FM,  noise-modulated  FM,  continuous  wave  (CW) 
interference,  linear-swept  CW,  log-swept  CW,  and  a  versatile  frequency-agile  pulsed- 
jamming  facility.  A  general  expression  for  the  complex  envelope  of  the  jamming 
waveform  is  given  in  terms  of  inphase  Ic(t)  and  quadrative  IJ(t)  components: 

I  (t)  =  Ic  ( t)-jlg  (t)  is  developed  as: 


(AW  +  W  ) 

Ic{t)  =  A[  1  +  \f^biCos{2TtfAMit)\  •  cos[ - -T--,  °tep  t 

+  ■7^E(ajsin(27r/FAf>0)  +  0)r 
JM  i 


(39a) 


and 


(A  W  +  W  ) 

lg  (f)  =  A[  1  4-  A^&^os^Tr/^f)]  •  sin[ - r—.-ep  t 


(39b) 


+  +  *]• 


A 

-  total  normalized  gain 

M 

-  modulation  index  of  AM  signal 

K 

-  weighting  term 

f  A  Mi 

-  AM  tone  frequency 

A  W 

-  offset  from  carrier  frequency 

w8lep 

-  sweep  step 

l(o.i) 

-  on/off,  logarithmic/linear  sweep  (I/O) 

A  f/fm 

-  modulation  index  FM  wave 
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weighting  term 
FM  tone  frequency 
phase  offset 


f  V 


'  v  - 


The  desired  jamming  scenarios  can  be  derived  from  the  above  equations  by 
substituting  zeros  for  the  unrelated  variables. 

For  the  simple  case  of  CW  interference,  the  following  variables  can  be  substituted 
with  the  constant  zero,  M,  6f,  fjMv  L,  // /m,  Aj,  and  fpyj  .  This  leaves: 

Ic(t)  =  AcOs(£Wt  +  9 ),  (40a) 


and 


/,  (t)  —  Asin(AWt  +  9),  (40b) 

which  are  the  formulas  for  CW  interference.  For  swept  CW  terms,  Wttep  and  tj0 
remain  in  the  original  equations  to  produce: 

IAW+  Wgt\ 

Ic(t)  =  Acos[ - T - ; - —  t  +  9]  (41a) 

C«(M) 

and 

t(AlV+  W.) 

Is  (0  =  Asin - : — 7 — —  t  +  9}  .  (4ib) 

C*(M) 

For  amplitude-modulated  interference,  the  following  variables  are  substituted  with 
zeros:  W 0tep  i^0  |j  ///m,  Aj,  and  Fp^j.  This  breaks  down  (39)  into  the  AM  standard 
form: 
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h  (0  =  A\\  +  M^,bt<x>s{2irfAMit)]cos{A  Wt  +  6)  (42a) 

i 

and 

/,(<)  =  A\\  +  M£blcos(2KfAMit)]sin{A  Wt  +  0)  .  (42b) 

t 

For  the  case  of  AM  noise,  the  term  btcos(2nf ^{t)  is  substituted  with  a 

I 

Gaussian  random  variable. 

FM  interference  can  be  simulated  when  M,  bit  ^»tepi  an<l  t(o,l)  are 

substituted  with  zeros.  The  equations: 

Ic  (t)  =  Acos[A  Wt  +  -^Ea;sin(27r/W)  +  °\  («a) 

Jm  j 

and 

h  (0  =  >lsin[A  Wt  +  ap'm{2KfFMjt)  +  0]  (43b) 

Jm  j 


result. 


A/ 

For  noise-modulated  FM,  the  quantity  — s — J]sin(27r/^yf)  -|-  0\  is  replaced  by  a 

Jm  j 

Gaussian  process. 


The  pulsed  frequency-agile  jammer  is  constructed  by  replacing  A  in  (39a  and  b) 


(44) 


A(<)  =  E  (t  -  tkyu_Y(t  -  tb  -  tk)e  jW>  . 

t=o 


Here  N{t)  is  the  same  counting  process  employed  in  the  impulsive  noise  channel 
with  an  additional  provision  for  perfectly  periodic  arrivals.  The  quantity  Wk  represents 
a  random  frequency  for  the  kt h  jamming  burst  and  is  drawn  from  a  distribution  under 
user  control.  The  random  jamming  amplitude  ak  is  chosen  from  either  a  power- 
Rayleigh  or  uniform  distribution. 

Although  only  these  five  distinct  choices  of  jamming  waveforms  are  available  as 
generic  jamming  options,  the  basic  jamming  module  allows  any  combination  of  the 
above  modulations  to  be  employed. 

Demodulator/  Receiver 

The  main  purpose  of  the  demodulator/receiver  function  module  is  to  perform  data 
demodulation  of  the  channel  output  r  ( t ).  That  is,  this  function  module  produces  the 
sequence  {.£,}  indicated  in  Figure  3-2,  which  is  applied  as  input  to  the  channel  decoder. 
The  demodulator/receiver  produces  an  output  once  per  baud  interval.  The  term  {£,} 
has  been  purposefully  allowed  to  be  a  vector  sequence  to  accommodate  various  soft- 
decision  decoding  schemes.  For  example,  in  an  M-ary  system,  the  receiver  might  well 
produce  an  ordered  list  of  the  /  largest  outputs  (/<A/)  during  successive  channel 
signaling  intervals.  This  is  the  so-called  list-of-f  quantization  scheme  [9)  which  provides 
a  reasonably  efficient  and  computationally  effective  means  for  implementing  soft- 
decision  decoding  in  a  variety  of  M-ary  signaling  situations. 

The  demodulator/receiver  module  contains  a  predetection  front  end  which 
precedes  the  data  demodulator  and/or  other  portions  of  this  module.  This  predetection 
section,  illustrated  in  Figure  3-12,  consists  of  a  zero-memory  nonlinear  (ZMNL)  device 
sandwiched  between  two  narrowband  filters.  The  filters  are  implemented  as  described 
previously  for  the  output  filtering  provided  in  the  modulator/transmitted  module.  A 
ZMNL  processing  capability  has  been  provided  in  the  ICS  to  allow  simulation  of 
intentional  nonlinear  processing,  such  as  might  be  used  to  mitigate  against  the  effects  of 


* 
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impulse  noise. 


The  various  ZMNL  processing  schemes  are  generally  of  the  limiting  variety  and 
have  been  implemented  in  the  ICS  as  illustrated  in  Figure  3-13.  Here  the  envelope 
/?0  (0  =  of  the  complex  envelope  s0  (f)  is  formed  and  passed 

through  an  instantaneous  or  zero-memory  nonlinear  (ZMNL)  device  with  input/output 
characteristic  </[/?0  (£)]  depending  only  upon  the  instantaneous  envelope  R0  (<). 
Forming  the  product 


«b  (0  =  rfflo  (Ofo  (0  («) 

allows  synthesis  of  a  variety  of  nonlinearities  of  the  limiting  or  saturation  type  by 
appropriate  choice  of  the  function  g  (R).  These  nonlinearities  are  all  characterized  by 
the  fact  that  they  depend  upon  the  envelope  and  not  upon  the  inphase  and  quadrature 
components  individually.  For  example,  the  choice  g  (R)  —  l/R  can  easily  be  shown 
to  correspond  to  an  ideal  bandpass  limiter.  The  choice 


fl  ;  0<R<1 

S(R)=\l/R  ;  /?>  1  ’  <47> 


on  the  other  hand,  represents  a  soft  limiter  or  clipper  characteristic.  Butterworth-type 
nonlinear  characteristics  of  various  cutoffs  and  orders  are  provided  under  user  control  in 
an  attempt  to  provide  a  realistic  and  yet  computationally  efficient  model  for 
demodulator/receiver  nonlinear  effects. 

It  should  be  noted  that  in  order  to  perform  its  basic  data  demodulation  functions, 
the  demodulator/receiver  module  must  perform  a  variety  of  ancillary  functions.  These 
ancillary  functions  include  symbol  synchronization  and  phase  tracking,  which  have  been 
implemented  for  BPSK  and  DPSK  modulation  strategies.  In  addition,  data  symbol  de¬ 
interleaving  is  performed  within  the  demodulator/receiver  module.  It  is  these  ancillary 
functions  which  seriously  complicate  the  design  of  this  function  module.  The  ICS  also 
has  the  capability  of  patching  symbol  synchronization  and  carrier  phase  and/or 
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frequency  from  transmitter  to  receiver.  This  facility  is  useful,  for  example,  when  one 
would  like  to  isolate  and  eliminate  the  various  synchronization  loop  tracking  dynamics 
from  consideration. 

Channel  Decoder 

The  channel  decoder  will  accept  the  sequence  {rj}  appearing  at  the  output  of  the 
demodulator/receiver  and,  through  appropriate  decoding  algorithms,  will  produce  the 
output  sequence  {fl,}  delivered  to  the  remote  destination.  Here  several  decoding  options 
are  provided  depending  upon  whether  block  or  convolutional  codes  are  used.  In  the 
latter  case,  the  Viterbi  Algorithm  (VA)  (cf.  [14],  [15])  is  employed  for  decoding  arbitrary 
convolutional  codes  with  constraint  lengths  in  the  range  3  <  K  <  9.  For  block  codes, 
several  hard-decisions  decoders,  using  essentially  table-look-up  techniques,  have  been 
implemented.  These  have  included  the  Hamming  (8,4)  and  the  Golay  (24,12)  codes 
which  are  typical  choices  in  a  wide  variety  of  applications.  In  addition,  a  fairly  general 
soft-decision  decoding  capability  for  arbitrary  block  codes  was  provided  using  the  Wolf 
algorithm  [16],  which  resembles,  in  many  ways,  the  VA.  Use  of  this  decoding  capability 
requires  explicit  knowledge  of  the  code-parity  check  matrix.  For  selected  codes  this 
information  is  stored  in  tables  on  the  host  computer.  The  user  then  specifies  a  code 
(e.g.,  the  Golay  (24,12)  code),  and  the  appropriate  parity-check  matrix  is  loaded  with 
the  decoder  program. 

3.3  OPERATION  OF  THE  ICS 

The  overall  technical  design  of  the  ICS  has  been  described  in  some  detail  in  the 
preceding  section.  In  this  section,  the  interactive  use  of  this  system  is  described. 
Careful  consideration  has  been  given  to  this  issue  to  allow  flexible  and  efficient  usage. 
The  general  philosophy  is  that  simulation  should  be  used  as  an  adjunct  to,  and  not  as  a 
substitute  for,  analysis.  The  communications-simulator  development  has  been 
consistent  with  this  philosophy.  While  the  basic  approach  is  Monte  Carlo  simulation, 
every  effort  has  been  made  to  incorporate  analytical  and  graphical  techniques  where  and 
when  appropriate  to  enhance  and  extend  the  capabilities  of  this  system. 


The  ICS  has  been  designed  to  be  used  in  four  distinct  modes  of  operation.  In  the 
first,  called  the  DESIGN  mode,  the  user  assembles  a  simulation  model  in  building-block 
fashion  from  basic  function  modules  provided  as  part  of  the  system.  Appropriate  user 
protocols  have  been  provided  to  guide  even  the  casual  user  through  the  DESIGN  stage. 
Specifically,  the  block  diagram  of  a  general  communications  system  (illustrated  in  Figure 
3-2)  is  displayed  on  the  graphics  CRT.  Using  the  cross-hair  curser,  the  user  selects 
successive  blocks,  and  a  menu  list  appears  describing  each  of  the  available  options  for 
each  functional  element.  The  user  chooses  from  this  list  by  entering  the  appropriate 
keyboard  response,  eventually  creating  a  DESIGN  file  representing  an  executable 
simulation  model  of  a  complete  end-to-end  communications  system.  This  DESIGN  file 
is  suitably  named  and  entered  into  the  system  library.  At  this  point,  the  user  can  either 
exit  from  the  DESIGN  mode  or  continue  creating  alternate  executable  DESIGN  files. 
This  is  the  most  straightforward  of  the  four  usage  modes  and  requires  no  further 
elaboration. 

The  second  usage  mode  is  called  the  VALIDATION  mode.  Here  the  user  is 
allowed  to  interactively  exercise  a  simulation  model  which  has  previously  been 
developed  in  a  DESIGN  session.  The  purpose  of  the  VALIDATION  mode  is  to  allow 
the  user  to  verify  a  simulation  model  before  a  commitment  to  time-consuming  Monte 
Carlo  simulation.  More  specifically,  the  user  is  allowed  to  view  time  waveforms  and/or 
frequency-domain  quantities  such  as  power  spectral  densities,  system  transfer  functions, 
etc.,  at  critical  points  in  the  system  to  be  simulated.  This  should  enable  the  user  to 
verify  the  end-to-end  behavior  of  the  simulation  model  for  the  communication  system 
under  study  or  the  input/output  behavior  of  any  specific  sub-module.  Situated  at  the 
graphics  terminal,  the  user  can  by  cross-hair  keyboard  entry  view  time-domain  or 
frequency-domain  quantities  as  if  a  scopeprobe  was  connected  to  the  corresponding 
point  in  an  actual  communications  system.  This  ability  is  termed  the  SCOPEPROBE 
concept.  This  capability  of  the  communication  simulator  should  not  be  considered  as 
merely  a  procedural  verification  stage  preceding  an  actual  simulation.  Rather,  the 
VALIDATION  mode  can  itself  be  employed  as  a  highly  flexible  and  useful  tool  in 
conceptual  communication-system  studies.  That  is,  by  employment  of  the 
VALIDATION  mode,  a  useful  and  realistic  appreciation  of  actual  system  behavior  can 
be  developed.  Such  an  appreciation  can  often  serve  as  a  useful  adjunct  to  detailed 


analysis.  For  example,  the  effects  of  various  filtering  operations  can  be  observed 
directly,  degradations  due  to  intersymbol  interference  or  multipath  can  be  assessed 
visually,  etc.  Some  typical  graphical  outputs  are  provided  in  the  next  subsection. 

In  order  to  modify  or  otherwise  change  a  simulation  model,  the  user  must  leave 
the  VALIDATION  stage  and  re-enter  the  DESIGN  mode.  Appropriate  editing 
capabilities  have  been  provided  to  allow  flexible  modification  of  previously  developed 
simulation  models.  At  any  rate,  the  user  is  at  some  point  ready  to  commit  a  simulation 
model  to  extensive  Monte  Carto  simulation.  It  is  at  this  point  that  the  user  enters  the 
third  stage,  called  the  SIMULATION  mode.  In  this  mode,  an  executable  DESIGN  file 
corresponding  to  a  simulation  model  must  be  identified  together  with  appropriate 
simulation  parameters.  Primary  interest  is  in  the  evaluation  of  bit-error  probability  Pb 
as  a  function  of  the  quantity  Ej N0.  Here  Eb  represents  the  signal  energy  per  bit  while 
Aq/2  represents  the  double-sided  noise  spectral  density  in  watts/Hz.  As  a  result,  the 
important  parameters  to  be  specified  are  both  an  initial  and  a  final  value  of  Eb/ N0  for 
which  simulation  results  are  desired  as  well  as  an  increment  in  Eb/N0.  In  addition, 
parameters  such  as  the  number  of  errors  to  be  collected  for  each  value  and  the 
maximum  number  of  iterations  must  be  specified.  These  latter  quantities  are  necessary 
in  order  to  ensure  statistical  validity  of  the  results.  The  communications  simulation  will 
then  loop  through  the  specified  range  of  Ej NQ.  This  is  the  only  parameter  for  which 
looped  operation  is  provided  in  the  SIMULATION  mode.  For  other  parameters  it  is 
required  that  an  appropriate  DESIGN  file  representing  a  simulation  model  with  specific 
parameter  choices  has  been  identified  for  subsequent  execution. 

In  the  SIMULATION  mode,  appropriate  data-logging  and  bookkeeping  software  is 
provided  to  allow  tabulation  of  bit-error  histories  and  evaluation  of  bit-error 
probabilities  as  a  function  of  Eb/N0.  This  software  is  considered  an  important  part  of 
the  communications  simulator. 

The  last  and  final  mode  of  operation  is  called  the  ANALYSIS  mode.  This  mode 
provides  two  basic  functions.  The  primary  function  is  to  provide  the  graphical  display 
interpretation  of  results  obtained  during  a  typical  SIMULATION  session.  Appropriate 
display  drivers  and  associated  software  have  been  provided  to  allow  results  to  be 


displayed  in  a  format  which  allows  easy  interpretation  by  communication  systems 
engineers.  For  example,  a  basic  requirement  is  to  display  curves  representing  simulated 
bit-error  probability  Pj  as  a  function  of  Eb/ Nq.  In  those  cases  where  closed-form 
analytical  expressions  or  bounds  on  Pb  vs.  EJNq  are  available,  these  are  included  as 
user-selected  options.  The  user  protocol  in  the  ANALYSIS  mode  will  provide  guidance 
to  the  user  concerning  any  such  facilities  available  for  the  system  configuration  and/or 
parameter  choices  under  study.  For  example,  fairly  tight  upper  bounds  are  available 
[14],  [17]  on  Pb  vs.  Ebf  Nq  for  convolutional  codes  in  conjunction  with  Viterbi  decoding 
and  employing  a  variety  of  modulation  formats  on  the  AWGN  channel.  Similar  bounds 
have  been  derived  [18]  for  soft-decision  decoding  of  block  codes  employing  the  Wolf 
algorithm.  These  bounds,  for  the  most  part,  neglect  the  effects  of  narrowband  filtering, 
imperfect  phase  and/or  symbol  synchronization,  intersymbol  interference,  etc. 
Nevertheless,  the  readily  computed  upper  bounds  are  useful  in  providing  a  context  or 
perspective  in  which  to  assess  the  degradations  due  to  these  effects.  For  this  reason,  the 
ability  to  compute  and  display  these  idealized  performance  bounds  together  with  results 
obtained  by  simulation  has  been  incorporated  as  a  useful  feature  of  the  ICS. 

In  the  course  of  developing  the  ICS,  it  was  found  to  be  useful  to  access  the 
computational  and  display  capabilities  associated  with  the  ANALYSIS  mode  on  a 
stand-alone  basis.  For  example,  the  ability  to  compute  and  plot  upper  bounds,  and  in 
some  cases,  actual  theoretical  performance  of  selected  communication  systems,  has 
proven  useful  in  a  number  of  related  system  studies.  As  a  result,  the  ANALYSIS  mode 
has  been  designed  to  allow  dual  entry,  thus  providing  a  secondary  function  as  a  stand¬ 
alone  computational  and  display  resource.  The  library  associated  with  this  mode  can  be 
easily  updated  to  provide  an  expanding  analysis  tool  quite  independent  of  the  simulation 
capabilities  associated  with  the  overall  system. 

A  block  diagram  illustrating  the  overall  software  organization  of  the 
communications  simulator  is  presented  in  Figure  3-14. 
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3.4  TYPICAL  GRAPHICAL  OUTPUT 


In  this  section,  some  of  the  hard-copy  graphical  outputs  obtained  with  the  ICS  are 
discussed.  First  considered  are  the  time-domain  waveforms  obtainable  in  the 
VALIDATION  mode.  Frequency-domain  graphical  display  and/or  output  capabilities 
are  available  and  will  be  described  below.  Virtually  any  time  waveform  can  be  observed 
individually  or  in  combination  with  other  related  waveforms.  The  user  must  first 
specify  a  time  or  viewing  window  consisting  of  a  fixed  number  of  baud.  A  segment  of  a 
waveform  is  then  displayed  within  the  viewing  window,  and  the  contents  can  be  indexed 
through  a  waveform  file  either  in  fixed  increments  or  continuously  scrolled  through  a 
large  file.  The  scrolling  can  be  stopped  and  restarted  under  keyboard  control. 

It  was  found  convenient  to  classify  time-domain  waveforms  into  five  broad, 
categories: 

1.  BASEBAND  Waveforms 

2.  CHANNEL  Waveforms 

3.  Eye  Diagrams 

4.  Continuous  I/Q  Plots 

5.  Sampled  I/Q  Plots. 

The  BASEBAND  Waveforms  include  all  time-domain  quantities  which  allow  Iowpass 
signal  representations  and  are  in  some  sense  synchronous  with  the  baud  rate.  This 
would  include  source  outputs,  source/channel  encoder  and  decoder  outputs,  and 
matched  filter  outputs  sampled  at  the  baud  rate.  By  comparison,  CHANNEL 
waveforms  include  all  time-domain  waveforms  which  require  complex  representations, 
either  in  polar  coordinates  or  in  terms  of  I/Q  components.  Actually,  the  user  protocol 
allows  the  user  to  specify  either  representation.  Frequency-domain  plotting  and  display 
capabilities  are  available  for  any  CHANNEL  waveforms  appearing  within  the  viewing 


window  and  using  a  wide  variety  of  possible  window  functions.  Detailed  information  on 
the  various  window  options  is  provided  in  the  user  protocol.  Again,  the  frequency- 
domain  representation  can  be  displayed  in  either  polar  or  rectangular  (i.e.  I/Q)  format 
under  user  control.  The  eye  diagram  represents  the  continuous  running  sum  of  the 
“undumped”  matched  filter  output.  The  continuous  I/Q  plot  is  a  polar  representation 
of  sampled  complex  data  which  is  available  at  the  output  of  any  module  generating 
complex  data.  Finally,  the  sampled  I/Q  plot  is  a  polar  representation  of  the  I  versus  Q 
matched  filter  output  statistic. 

Some  typical  BASEBAND  signaling  waveforms  are  illustrated  in  Figure  3-15  for  a 
K=6,  R=l/2  eonvolutionally  encoded  QPSK  system  on  the  AWGN  channel  with 
EJNq  =  2 OdB.  In  this  case  the  observation  window  consists  of  8  baud,  and  a 
MULTIPLE  WAYUFORM  plotting  option  is  used.  More  specifically,  while  the  user  can 
look  at  any  individual  BASEBAND  waveform  in  the  COMPOSE  plotting  and  display 
mode,  the  ICS  provides  the  capability  of  plotting  certain  logical  combinations  in  the 
MULTIPLE  WAVEFORM  plotting  option.  These  figures  illustrate  the  source  output 
and  Viterbi  decoder  output  in  time  alignment  in  the  right-hand  column,  while  the 
encoder  output  and  matched  filter  output  are  in  time  alignment  along  the  left-hand 
column.  This  is  a  logical  arrangement  of  multiple  BASEBAND  waveforms  and  is 
provided  to  the  user  as  a  display  option. 

Typical  CHANNEL  signaling  waveforms  are  illustrated  in  Figure  3-16.  The 
waveforms  in  Figure  3-16  correspond  to  a  coherent  QPSK  system  operating  on  the 
AWGN  channel  at  the  relatively  large  value  EJNq  =  2 OdB.  Here  again,  use  is  made 
of  a  MULTIPLE  WAVEFORM  plotting  capability  with  the  left-hand  column  associated 
with  the  modulator/transmitter  output  s  ( t ),  while  the  right-hand  column  is  associated 
with  the  channel  output  r  (<).  Furthermore,  the  CHANNEL  waveforms  in  this  case  are 
plotted  in  polar  coordinates  in  terms  of  envelope  and  phase.  Similar  plots  in  rectangular 
coordinates  using  I/Q  components  could  have  been  chosen  under  user  control.  Likewise, 
individual  CHANNEL  waveforms  could  have  been  selected  under  the  COMPOSE 
plotting  and  display  option.  In  this  case,  frequency-domain  representations  could  be 
computed  and  displayed. 
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For  a  complete  description  of  the  remaining  graphical  display  options,  see  the  ICS 
User’s  Manual. 

A  number  of  additional  graphical  plotting  and  display  capabilities  exist  in  the 
ANALYSIS  mode.  For  the  most  part,  these  consist  of  accessing,  formatting,  and 
displaying  error  statistics  or  noise  amplitude  probability  distributions  (APDs)  generated 
in  executing  an  appropriate  DESIGN  file  in  the  SIMULATION  mode.  For  example,  a 
DESIGN  file,  representing  a  coherent  QPSK  system  operating  on  the  AWGN  channel 
and  employing  a  K=6,  R=l/2  convolutional  code,  was  executed  in  the  SIMULATION 
mode,  and  both  raw  and  decoded  error  statistics  were  collected.  Figure  3-17  illustrates 
the  measured  symbol-error  probability  P€  and  the  decoded  bit-error  probability  Pb. 
The  decoded  bit-error  probability  is  plotted  as  a  function  of  EJ Nq,  while  the  symbol- 
error  probability  is  plotted  as  a  function  of  Ej Nq.  Here,  Et  is  the  signal  energy  on  a 
per-baud  basis  and  is  related  to  according  to  Et  =  RE j,  where  R  is  the  code  rate 
(in  this  case  R=l/2).  Also  included  in  the  respective  plots  in  Figure  3-17  are  the 
theoretical  symbol-error  probability  Pe  and  the  computed  upper  bound  on  Py  The 
simulation  results  compare  favorably  to  these  computed  results.  These  results  provide  a 
complete  picture  of  the  performance  of  the  coded  coherent  QPSK  system  on  the  AWGN 
channel.  The  total  elapsed  time  to  execute  the  SIMULATION  mode  (in  this  case  over 
the  range  0 dB  <  EJ N0  <  6 dB  in  steps  of  0.5  dB)  was  2  hours.  This  is  significant 
when  one  observes  that  error  events  as  rare  as  one  in  107  have  been  accurately 
measured.  Comparable  simulation  capabilities  would  be  impossible  in  existing  general- 
purpose  machines.  Indeed,  it  is  questionable  whether  such  an  extensive  simulation 
would  be  undertaken  in  the  first  place. 

Other  graphical  output  available  in  this  case  would  include  histograms  of  the 
additive  channel  noise  amplitude  and  phase,  two-dimensional  histograms  of  the  phase 
and  time  synchronization  errors,  and  three-dimensional  plots  of  the  fading  Channel 
scattering  function.  This  scattering  function  is  shown  in  Figure  3-18  for  an  21-tap  delay 
line  model. 
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4.  CONCL USIONS  AND  RECOMMENDA  TIONS 


The  enhanced  ICS  provides  a  comprehensive  facility  for  the  modeling,  simulation, 
and  analysis  of  point-to-point  digital  communications  systems.  Applications  modules 
available  in  ICS  include  modulators/demodulators,  coders/decoders,  filters, 
nonlinearities,  and  a  comprehensive  set  of  modules  for  the  simulation  of  realistic 
transmission  impairments.  These  include  fading  dispersive  phenomena  as  well  as  ECM 
interference.  The  ICS  is  being  employed  in  the  evaluation  of  single-link  communication 
systems  and  in  the  development  of  state-of-the-art  signal  processing  concepts  for 
military,  space,  and  commercial  applications. 

The  major  difficulties  encountered  in  the  ICS  ENHANCEMENTS  efforts  were 
caused  by  hardware/software  limitations  of  the  target  computer  configuration,  a  PDP- 
11/40  computer  residing  at  DICEF.  ICS  has  become  a  large-scale  software  system  and, 
as  such,  will  run  efficiently  only  on  a  virtual  memory  machine  such  as  a  VAX-11  series 
computer.  The  limited  task  size  (32K)  of  PDP-11  series  computers  severely  complicates 
program  development  by  requiring  complex  overlay  structures.  Since  overlays  run  much 
less  efficiently  than  single  image  tasks,  the  overall  speed  and,  hence,  usefulness  of  a 
PDP-11  based  ICS  is  also  severely  compromised.  The  replacement  of  DICEF  PDP- 
11/40  computer  with  a  VAX-11  series  equivalent  computer  is  recommended  as  the 
highest  priority  for  future  system  improvements. 

The  applications  modules  provided  in  the  current  enhancement  effort  addressed 
channel  modeling,  spread-spectrum  techniques,  ECM  waveform  generation,  and 
enhanced  graphics  capabilities.  The  major  enhancement  which  was  planned  but  not 
implemented  was  an  adaptive  array  simulation  capability.  The  usefulness  of  the  ICS 
would  be  greatly  enhanced  by  the  inclusion  of  such  a  capability.  Other  improvements 
which  would  be  useful  include  improved  tutorial  interaction,  adaptive  equalizers,  and 
more  comprehensive  block  and  concatenated  coding  schemes.  The  capability  to  collect 
bit-error  and  symbol-error  statistics  other  than  average  error  probability  would  prove 
quite  useful  on  fading  and/or  burst  noise  channels.  Any  of  these 
applications/enhancements  should  be  accomplished  only  after  an  upgrade  to  a  virtual 
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APPENDIX  A 


BLOCK  CODES 


In  a  block  encoder,  the  binary  sequence  {aj}  is  segmented  into  blocks  of  size  k  to 
which  n-k  redundant  bits  are  added  to  produce  a  rate  R=k/n  code  measured  in  units  of 
information  bits  per  transmitted  channel  symbol.  There  is  complete  independence 
between  blocks.  The  set  of  all  binary  n-tuples  is  a  vector  space.  A  set  of  these  vectors 
of  length  n  is  called  a  linear  block  code  if  and  only  if  it  is  a  subspace  of  the  vector  space 
of  n-tuples.  In  the  binary  case,  and  in  fact  over  a  field  of  p  elements  where  p  is  prime, 
every  group  of  vectors  is  also  a  subspace.  The  term  group  code  is  common  terminology 
for  binary  linear  codes. 

The  Hamming  weight  of  a  vector  X,  denoted  tfl(  j)  is  defined  to  be  the  number  of 
nonzero  components.  Since  the  Hamming  distance  between  two  vectors  X^  and  Xq  is  the 
number  of  positions  in  which  they  differ,  the  distance  between  x^  and  X_2  is  equal  to 
w(  -  x2y  If  Xx  and  X2  are  both  code  words  of  a  linear  block  code,  xx  -  Xq  must  also 
be  a  code  word,  since  the  set  of  all  code  words  is  a  vector  space.  Therefore,  the  distance 
between  any  two  code  vectors  equals  the  weight  of  some  other  code  vector,  and  the 
minimum  distance  for  a  linear  code  equals  the  minimum  weight  of  its  nonzero  vectors. 
This  property  is  extremely  helpful  in  analyzing  the  error-correction  capabilities  of  linear 
codes. 


Nearly  all  known  block  and  convolutional  codes  are  based  on  the  concept  of  Ham¬ 
ming  distance.  Any  set  of  basis  vectors  for  a  linear  block  code  X  can  be  considered  as 
rows  of  a  matrix  G,  called  a  generator  matrix  of  X.  The  row  space  of  G  is  the  linear 
code  X,  and  the  vector  is  a  code  vector  if  and  only  if  it  is  a  linear  combination  of  the 
rows  of  G.  If  the  dimension  of  the  vector  space  X  is  k,  the  number  of  rows  of  G 
(which  equals  the  rank  of  G,  since  the  rows  must  be  linearly  independent)  is  k.  If  any 
two  linear  combinations  were  equal,  there  would  be  a  dependence  relation  among  rows 
of  G.  Thus,  each  distinct  linear  combination  gives  a  distinct  code  vector,  and,  since 


there  are  k  coefficients  with  2  possible  values  for  each,  there  are  2  code  vectors  in  X. 
Such  a  code  is  called  an  (n,k)  code. 


Unless  k  is  small,  the  matrix  description  is  much  more  compact  than  a  list  of  code 

g 

vectors.  A  (50,  30)  code  is  described  by  a  30  x  50  matrix  but  has  more  than  10  code 
vectors. 


Example:  A  typical  code  with  N=7  and  K=4  has  a  generator  matrix  given  by 


G  = 


10  0  0 
0  10  0 
0  0  10 
0  0  0  1 


0  1  1 
1  0  1 
1  1  0 
1  1  1 


This  is  the  Hamming  (7,  4)  Code  with  minimum  distance  d=3. 


There  is  an  alternative  description  by  matrices.  Again,  if  X  is  a  sub-space  of  dimension 
k,  its  null  space  is  a  vector  space  X  of  dimension  n-k.  A  matrix  H  of  rank  n  -  k  whose 
row  space  is  X  can  be  made  with  a  basis  for  X  as  rows.  Then  X  is  the  null  space  of  X  , 
and  a  vector  x  is  in  X  only  if  it  is  orthogonal  to  every  row  of  H,  that  is,  if  and  only  if 


(A-l) 


If  X  =  (Xj,  Xg,  X^,  .  .  .,  Xn)  and  the  element  in  the  ith  row  and  jth  column  of  H 
is  denoted  h—.  Equation  A-l  implies  that  for  each  i  (that  is  for  each  row  of  H), 

E  =  xihij  =  0  •  (A*2) 

i 

Thus,  the  meaning  of  Equation  A-l  is  that  the  components  of  X  must  satisfy  a  set 
of  n  -  k  independent  equations.  Of  course,  any  linear  combination  of  the  Equation  A-2 
also  gives  an  equation  that  the  components  of  X  must  satisfy,  and  this  corresponds  to 


the  fact  that  X  is  orthogonal  to  every  vector  of  X .  These  equations  are  called 
generalized  parity  checks,  since  in  the  binary  case  they  are  simply  checks  for  even  parity 
on  certain  sets  of  symbols  in  the  code  word.  That  is,  for  each  row  of  H,  the  number  of 
l’s  in  X  that  corresponds  to  l’s  in  that  row  of  H  is  even  (for  the  binary  case)  if  and 
only  if  X  satisfies  Equation  A-l.  The  matrix  H  is  called  a  parity-check  matrix  of  X. 

Equation  A-l  holds  for  every  vector  x  in  X.  In  particular,  it  holds  for  the  k  basis 
vectors  of  the  matrix  G  These  k  equations  can  be  expressed  as  follows,  with  0 
denoting  the  kx  (n-k)  all-zero  matrix: 

f 

Example:  The  null  space  X  of  the  vector  space  X  of  the  previous  example  consists  of 
all  possible  linear  combinations  of  the  rows  of 

0  11110  0 
H  =  10  110  10. 

110  10  0  1 

♦ 

The  code  X  is  the  null  space  of  this  matrix.  For  each  vector  in  X  there  is  an  equation 
that  the  components  of  every  code  vector  must  satisfy.  For  example,  corresponding  to 
(OllllOl)is  the  equation 

OXj  ■+■  IA2  +  1X3  •+•  1X4  +  lXg  +  0X6  +  OA7  =  0 

that  must  be  satisfied  by  every  code  vector  (Xj,  X2,  X3,  X4,  X5,  Xg,  X7).  For  binary 
codes,  this  is  equivalent  to  having  an  even  number  of  l’s  among  the  last  four 
components  or  an  even  parity  check  on  the  last  four  components.  Note  that,  unlike 
vectors  over  the  field  of  real  numbers,  a  vector  can  be  orthogonal  to  itself.  For  example, 
(0  1  1  1  1  0  0)  is  in  both  X  and  X  . 


APPENDIX  B 


NARROWBAND  FILTERING 

In  implementing  a  digital  filter  on  a  digital  computer  or  with  special-purpose 
hardware,  the  input-output  relation  must  be  converted  to  a  computational  algorithm. 
The  algorithm  is  essentially  specified  in  terms  of  a  set  of  basic  computations  or  ele¬ 
ments.  For  the  implementation  of  discrete-time  systems  described  by  linear  constant- 
coefficient  difference  equations,  it  is  convenient  to  choose  as  these  elements  the  basic 
operations  of  addition,  delay,  and  multiplication  by  a  constant.  The  computational 
algorithm  for  implementing  the  filter  is  then  defined  by  a  structure  or  network  consist¬ 
ing  of  an  interconnection  of  these  basic  operations.  As  an  illustration,  consider  a  system 
with  a  system  function  of  the  form 


M 

t— 0 
N 

1-  V  V’  fc 

t=0 


ATz) 


The  difference  equation  relating  input  and  output  is  easily  written  down  directly  from 
the  system  function  and  is  given  by 

N  M 

y(n)  =  £  «*!/(«-*)  +  £  M 


We  can  interpret  Equation  B-2  directly  as  a  computational  algorithm  in  which  the 
delayed  values  of  the  input  are  multiplied  by  the  coefficients  bt~  ,  the  delayed  values  of 
the  output  are  multiplied  by  the  coefficients  a ,  and  all  the  resulting  products  are 


A  block-diagram  representation  for  the  general  difference  equation  of  Equation  B-2 
is  given  in  Figure  B-l.  The  network  of  Figure  B-l  is  an  explicit  graphical  representation 
of  the  difference  equation  B-2. 

In  general,  a  linear  shift-invariant  system  may  have  a  unit-sample  response  that  is 
of  finite  duration,  br  it  may  be  of  infinite  duration.  Because  of  the  properties  of  some 
digital  processing  techniques,  it  is  useful  to  distinguish  between  these  two  classes.  If  the 
unit-sample  response  is  of  finite  duration,  it  will  be  referred  to  as  a  finite  impulse 
response  (FIR)  system;  if  the  unit-sample  response  is  of  infinite  duration,  it  will  be 
referred  to  as  an  infinite  impulse  response  (HR)  system.  If  N  =  0  in  Equation  B-2,  so 
that 

M  =  ~r  E  MM  - 

“0  r=0 

it  corresponds  to  an  FIR  system.  In  fact,  the  above  difference  equation  is  identical  to 
the  convolution  sum,  and  hence  it  follows  directly  that 

X 

^  n~  0,1,..., A/ 

0,  otherwise 

An  FIR  system  can  always  be  described  by  a  difference  equation  of  the  form  of  Equation 
B-2  with  N  =  0.  In  contrast,  for  an  IIR  system,  N  must  be  greater  than  zero. 

Design  of  IIR  Filters 

The  traditional  approach  to  the  design  of  IIR  digital  filters  involves  the 
transformation  of  an  analog  filter  into  a  digital  filter  meeting  prescribed  specifications. 


B-2 


The  most  common  design  procedure  is  based  on  integrating  the  differential 
equation  and  then  using  a  numerical  approximation  to  the  integral.  This  procedure 
leads  to  the  so  called  bilinear  transformation  where  the  z  transforms  of  the  digital 
filter’s  impulse  response  H(z)  are  obtained  from  the  Laplace  transform  H(s)  of  the 
corresponding  analog  filter  by  the  substitution 


That  is, 


y^)  =  flW|s=(2/n(l-r')(l+2->). 


Solving  Equation  B-3  for  Z  gives: 

_  l+(T,/2)s 

l-(r,/2).s  '  (B-5' 

To  demonstrate  that  this  mapping  has  the  property  that  the  imaginary  axis  in  the 
s-plane  maps  onto  the  unit  circle,  consider  Z=eJW  .  Then,  from  Equation  B-5,  S  is 
given  by 


5=  JL  l-e~3W 

T  i+e-v 


T  cos (cj/s) 


Thus  for  z  on  the  unit  circle,  <7=0  and  Q  and  u)  are  related  by 

7—-  =  tan(a;/2)  . 

This  relationship  is  plotted  in  Figure  B-2.  From  the  figure  it  is  clear  that  the  positive 
and  negative  imaginary  axes  of  the  s-plane  are  mapped,  respectively,  into  the  upper  and 
lower  halves  of  the  unit  circle  in  the  z-plane. 

In  addition  to  the  fact  that  the  imaginary  axis  in  the  s-plane  maps  to  the  unit 
circle  in  the  z-plane,  the  left  half  of  the  s-plane  maps  to  the  inside  of  the  unit  circle  and 
the  right  half  of  the  s-plane  maps  to  the  outside  of  the  unit  circle,  as  depicted  in  Figure 
B-3. 

This  can  be  seen  by  referring  to  Equation  B-5:  for  the  real  part  of  s  negative,  the 
magnitude  of  the  factor  (l+sT/2)  /  ( \-sT/2 )  is  less  than  unity,  corresponding  to  the 
inside  of  the  unit  circle.  Conversely,  for  the  real  part  of  s  positive,  the  magnitude  of 
that  ratio  is  greater  than  unity,  corresponding  to  the  outside  of  the  unit  circle.  Thus  we 
see  that  use  of  the  bilinear  transformation  yields  stable  digital  filters  from  stable  analog 
filters. 

Digital  Butterworth  Filters 

Butterworth  filters  are  defined  by  the  property  that  the  magnitude  response  is 
maximally  flat  in  the  passband.  For  an  Nth-order  lowpass  filter,  this  means  that  the 
first  2N  -  1  derivatives  of  the  squared  magnitude  function  are  zero  at  fi  =  0  .  Another 
property  is  that  the  approximation  is  monotonic  in  the  passband  and  the  stopband. 
The  squared  magnitude  function  for  an  analog  Butterworth  filter  is  of  the  form 


1 


i+(jn/jnc)2N 


(B-6) 


As  the  parameter  N  in  Equation  B-6  increases,  the  filter  characteristics  become 
sharper;  that  is,  they  remain  closer  to  unity  over  more  of  the  passband  and  become  close 
to  zero  more  rapidly  in  the  stopband,  although  the  magnitude  function  at  the  cutoff 
frequency  f2c  will  always  be  l/\/2  because  of  the  nature  of  Equation  B-6.  The 
dependence  of  the  Butterworth  filter  characteristic  on  the  parameter  N  is  indicated  in 
Figure  B-4. 


From  the  squared  magnitude  function  in  Equation  B-6,  we  observe  that 
Ha(s)Ha(-s)  must  be  of  the  form 


HMHJl-s) 


1 _ 

i+(5/jnc)“  ' 


(B-7) 


The  roots  of  the  denominator  polynomial  (the  poles  of  the  squared  magnitude  function) 
are  then  at 


*,  =  (-»^v(jn.)  • 

Thus,  there  are  2N  poles  equally  spaced  in  angle  on  a  circle  of  radius  Oc  in  the  s-plane. 
The  poles  are  symmetrically  located  with  respect  to  the  imaginary  axis.  A  pole  never 
falls  on  the  imaginary  axis,  and  one  occurs  on  the  real  axis  for  n  odd  but  not  for  N  even. 
The  angular  spacing  between  the  poles  on  the  circle  is  7T / N  radians.  In  designing  a 
Butterworth  filter  using  the  bilinear  transformation,  the  most  straightforward  procedure 
is  to  first  determine  the  location  of  the  poles  in  the  s-plane  and  then  to  map  the  left- 
hand  plane  poles  to  the  z-plane. 


B-8 


Digital  Chebyshev  Filters 


In  a  Butterworth  filter,  the  frequency  characteristic  is  monotonic  in  both  the 
passband  and  the  stopband.  Consequently,  if  the  filter  specifications  are  in  terms  of,  let 
us  say,  maximum  passband  approximation  error,  the  specifications  are  exceeded  toward 
the  low-frequency  end  of  the  passband.  A  more  efficient  approach,  which  usually  leads 
to  a  lower-order  filter,  is  to  distribute  the  accuracy  of  the  approximation  uniformly  over 
the  passband  or  the  stopband  or  both.  This  is  accomplished  by  choosing  an 
approximation  that  has  an  equiripple  behavior  rather  than  a  monotonic  behavior.  The 
class  of  Chebyshev  filters  has  the  property  that  the  magnitude  of  the  frequency  response 
is  either  equiripple  in  the  passband  and  monotonic  in  the  stopband  or  monotonic  in  the 
passband  and  equiripple  in  the  stopband. 


The  analytic  form  of  the  squared  magnitude  function  for  the  Chebyshev  filter  is 


HM 


2 


1 

i+£2F2"(n/nc)  ’ 


(B-8) 


where  V x)  is  the  Nth-order  Chebyshev  polynomial  defined  as 
Vyv(x)  =  cos(Ncos~lx)  . 


(B-9) 


For  example,  for  N  =  0,  V jj(x)=l~  ;  for  iV=l ,  VyV(x)=COs(cos  1x)=x~  ;  for 
/V=2,  ViV(j)=cos(2cos“1x)=2x2-l~  ;  etc. 

From  Equation  B-9,  which  defines  the  Chebyshev  polynomials,  it  is 
straightforward  to  obtain  a  recurrence  formula  from  which  V}y+1(ar)  can  be  obtained 
from  Vf\,(x)  and  by  applying  trigonometric  identities  to  Equation  B-9,  with 

the  result  that 


VW*)  =  %*vfk*)~vN-\(x)  ■ 


(B-10) 


From  Equation  B-9  we  note  that  V2 pj(x)  varies  between  zero  and  unity  for  x  between 
zero  and  unity.  For  x  greater  than  unity,  COS-1  x  is  imaginary,  so  that  V jy(x)  behaves 
as  a  hyperbolic  cosine  and  consequently  increases  monotonically  for  x  greater  than 
unity.  Referring  to  Equation  B-8,  then,  ripples  between  1  and  l/(l-f£2)' 

for  0<=nnc<=  1  and  decreases  monotonically  for  f2/Oc>l.  Three  parameters  are 
required  to  specify  the  filter:  £,  flc  ,  and  N.  In  a  typical  design,  £  is  specified  by  the 
allowable  passband  ripple,  and  Qc  is  specified  by  the  desired  cutoff  frequency.  The 
order  N  is  then  chosen  so  that  the  stopband  specifications  are  met. 

The  poles  of  the  Chebyshev  filter  lie  on  an  ellipse  in  the  s-plane. 

Design  of  FIR  Filters 

The  previous  section  has  been  concerned  with  the  design  of  OR  filters  using  the 
bilinear  transformation.  Although  such  filters  have  many  attractive  features,  they  also 
have  a  number  of  disadvantages.  For  example,  if  one  wishes  to  take  advantage  of  the 
computational  speed  of  an  FFT  implementation,  a  finite-duration  impulse  response  is 
essential.  Also,  IIR  filters  generally  achieve  excellent  amplitude  response  at  the  expense 
of  non-linear  phase.  In  contrast,  FIR  filters  can  have  exactly  linear  phase.  Thus,  design 
techniques  for  FIR  filters  are  of  considerable  interest. 

The  system  function  of  a  causal  FIR  filter  is  of  the  form 

H{z)  =  *£  V)e->"  , 

n=0 

i.e.,  H  .(z)  is  a  polynomial  in  of  degree  N  -  1.  Thus  H  (z)  has  N  -  1  zeros  that  can  be 
located  anywhere  in  the  finite  z-plane  and  N  -  1  poles  all  of  which  lie  at  z  =  0.  The 


frequency  response  H(erfis  the  trigonometric  polynomial 


H(e^)=  N'£  h(n)e-»}n  . 

R=0 


(B-ll) 


We  recall  that  any  finite-dqration  sequence  is  completely  specified  by  N  samples  of  its 
Fourier  transform,  so  that  the  design  of  an  FIR  filter  may  be  accomplished  by  finding 
either  its  impulse  response  coefficients  or  N  samples  of  its  frequency  response. 

If  the  impulse  response  satisfies  the  condition 

h(n)  =  h(N-l-n)  ,  (B-12) 


the  filter  has  linear  phase.  This  is  easily  shown  by  substituting  Equation  B-12  into  B- 
11,  thereby  obtaining 


He**)  = 


JV-l  (N-3)/2 


-M(N- 1)/2) 


^/2  2 h(n)  cos  ju|  j ,  N  odd  (B-13A) 


N  even  (B-13B) 


It  is  seen  from  these  equations  that  the  condition  of  Equation  B-12  implies  a  linear 
phase  shift  corresponding  to  a  delay  of  (N  -  l)/2  samples.  We  note  that  for  the  case  of 
N  odd,  the  phase  shift  corresponds  to  an  integer  number  of  samples  delay,  while  for  N 
even,  the  delay  is  an  integer  plus  one-half  sample. 

Since  linear  phase  is  generally  desirable  and  often  a  necessity,  and  since  linear 
phase  often  simplifies  a  design  procedure,  our  discussion  will  center  on  linear-phase 


The  most  straightforward  approach  to  FIR  filter  design  is  to  obtain  a  finite-length 
impulse  response  by  truncating  an  infinite-duration  impulse  response  sequence.  If  we 
suppose  that  is  an  ideal  desired  frequency  response,  then 


Hde»=  £ 


(B-14A) 


where  hj(n)  is  the  corresponding  impulse  response  sequence,  i.e., 


(B-14B) 


In  general,  for  a  frequency  selective  filter  may  be  piecewise  constant  with 

discontinuities  at  the  boundaries  between  bands.  In  such  cases  the  sequence  hj[n)  is  of 
infinite  duration,  and  it  must  be  truncated  to  obtain  a  finite-duration  impulse  response. 
As  we  have  pointed  out  before,  B-14  can  be  thought  of  as  a  Fourier  series  representation 
of  the  periodic  frequency  response  e^),  with  the  sequence  hj^  n)  playing  the  role  of 
the  “Fourier  coefficients.”  Thus,  the  approximation  of  an  ideal  filter  specification  by 
truncation  of  the  ideal  impube  response  b  identical  to  the  study  of  the  convergence  of 
Fourier  series,  a  subject  that  has  received  a  great  deal  of  study  since  the  middle  of  the 
eighteenth  century.  The  most  familiar  concept  from  thb  theory  b  the  Gibbs 
phenomenon. 

If  hj^n)  has  infinite  duration,  one  way  to  obtain  a  finite-duration  causal  impube 
response  b  to  simply  truncate  h(n),  i.e.,  define 


1/  \  _  M«)>  0  <  =  n  <=  N- 1 
H\n)  —  o,  otherwise 


(IMS) 


In  general,  we  can  represent  h(n)  as  the  product  of  the  desired  impulse  response  and  a 
finite-duration  “window”  w(n),  i.e., 


h(n)=hj[n)ii^n)  . 


(B-10) 


By  tapering  the  window  smoothly  to  zero  at  each  end,  the  height  of  the  sidelobes 
in  H(e **')  can  be  decreased  at  the  expense  of  a  wider  main  lobe.  Typical  windows  for 
this  application  can  be  found  in  [B-2]. 

Another  popular  method  of  designing  FIR  filters  is  by  sampling  the  frequency 
response  of  the  desired  analog  filter.  A  finite  duration  sequence  can  be  represented  by 
its  discrete  Fourier  transform. 


Thus  an  FIR  filter  has  a  representation  in  terms  of  the  “frequency  samples” 


H  (k)  —  H(z) 


N-l 


n=0 


H(z)  can  be  represented  in  terms  of  the  samples  H(k)  by  the  expression 


H(z)  = 


1  -zN  y  H  (k) 

N  UeA2«iml  ■ 


(B-17) 


Equation  B-17  serves  as  the  basis  of  the  frequency-sampling  realization  of  an  FIR  filter. 
If  we  let  z=e}0J,  the  frequency  response  has  the  representation 


H[e^)  = 


1  -e~juN  N-l 

£0  1  -e’&'IWe-i" 


n  W 


N 


(B-18) 


_  e-M(N-i)/2)  N-i  ~  sin  [N[u  -  {2Tr/N)k)/2] 

N  sin  [(a;  -  (2n/N)k)/2\ 


k=0 


B-14 


Equation  B-18  suggests  a  simple  but  rather  naive  approach  to  filter  design,  i.e.,  to 
specify  the  filter  in  terms  of  samples  of  one  period  of  the  desired  frequency  response 

relying  on  the  interpolation  indicated  in  Equation  B-18  to  “fill  in  the  gaps”  in  the 
frequency  response. 

Frequency  sampling  designs  are  particularly  attractive  for  narrow-band  frequency 
selective  filters  where  only  a  few  of  the  samples  of  the  frequency  response  are  nonzero. 
In  such  cases,  a  frequency  sampling  realization  may  be  considerably  more  efficient  than 
either  direct  convolution  or  convolution  using  the  DFT.  In  general,  even  if  more  than  a 
few  samples  are  nonzero,  the  frequency-sampling  design  method  yields  excellent  results. 

Although  the  results  of  the  frequency  sampling  design  procedure  can  be  very  good, 
it  is  still  not  optimum  in  a  mini-max  sense.  We  also  lack  flexibility  in  specifying  the 
stopband  and  passband  frequencies  WB  or  Wp  ,  since  all  samples  are  constrained  to 
2ttK/N.  It  is  possible  to  design  optimal  (in  a  prescribed  sense)  FIR  digital  filters. 
These  filters  are  called  equiripple  because  of  the  form  of  their  magnitude  response.  If 
the  approximation  error  in  the  digital  filter’s  frequency  response  were  spread  out 
uniformly  in  frequency,  a  given  design  specification  could  be  met  with  a  lower-order 
filter  than  if  the  approximation  just  meets  the  specification  at  certain  frequencies  and 
far  exceeds  it  at  others.  Suppose  a  digital  filter  design  is  desired  in  accordance  with  the 
tolerance  scheme  of  Figure  B-5.  That  is,  it  is  desired  to  approximate  1  in  the  frequency 
range  W,<=|  W]<=  Wp  within  a  tolerance  6l—K62  and  it  is  desired  to 
approximate  zero  in  the  frequency  range  W,<  =  |  VF|  <=7T  .  The  parameter  K  controls 
the  magnitude  of  the  ripple  in  the  passband  relative  to  the  magnitude  of  the  ripple  in 
the  stopband.  Parks  and  McClellan  [B-3]  have  developed  an  iterative  procedure  to 
minimize  the  value  of  62  given  the  filter  order  fi,  the  value  of  the  stopband  frequency 
Wr  the  value  of  the  passband  cutoff  frequency  Wt  ,  and  the  weighting  ratio  K.  It  can 
be  shown  that  this  procedure  results  in  a  filter  with  a  frequency  response  which  exactly 


B>15 


meets  the  maximum  magnitude  error  specification  of  exactly  M  +  2  frequencies 
including  Wa  and  WP  .  This  filter  is  optimum  in  the  sense  that  no  other  linear  phase 
FIE  filter  of  the  same  order  can  have  a  smaller  value  of  maximum  error  in  the  passband 
and  stopband. 

The  Parks  -  McClellan  algorithm  is  incorporated  in  the  ICS  design  software.  A 
typical  response  for  the  case 

M=  64,  Wp=.  5tt,  W=£it 

and  K  =  1  is  illustrated  in  Figure  B-6.  In  this  case,  the  frequency  axis  is  scaled  relative 
to  the  baud  rate,  and  thus  the  frequency  response  is  plotted  from  ~Ng/2  to  NJ 2  . 
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